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Abstract 

We use the Berry-phase-based theory of macroscopic polarization of dielec- 
tric crystals formulated in terms of Wannier functions, and state-of-the-art 
Gaussian basis functions, to obtain benchmark ab initio Hartree-Fock values 
of the Born effective charges of ionic compounds LiH, LiF, LiCl, NaF, and 
NaCl. We find excellent agreement with the experimental values for all the 
compounds except LiCl and NaCl, for which the disagreement with the exper- 
iments is close to 10% and 16%, respectively. This may imply the importance 
of many-body effects in those systems. 
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I. INTRODUCTION 



The Born effective charge (also called transverse or dynamic charge) of a crystalline 
system, defined as the induced polarization due to a unit sublattice displacement!^ is a 
fundamental quantity connecting the electrostatic fields of the lattice, to its phononic prop- 
erties. It contains important information not only about the electronic structure and the 
bonding properties of the system, but also about the coupling of its longitudinal and trans- 
verse optical phonon modes to the external infrared radiationBi^ In addition, the Born 
charges also find use in first-principles-based construction of effective Hamiltonians describ- 
ing phase transitions in ferroelectric materials.^! Traditionally, the ab initio computations of 
Born charges of dielectrics have been performed either within the density-functional linear- 
response theory,i0 or the density-functional perturbation theory.i0 At a morephenomeno- 
logical level, the bond-orbital method of Harrison has been very insightfulca Recently, 
however, a very elegant formalism has been proposed by King-Smith and VanderbiltJI^I 
and Resta. which formulates the general problem of symmetry breaking induced by the 
macroscopic polarization (of which the Born charge is a special case) of a crystalline dielec- 
tric, in terms of the Berry phase of its wave function. This Berry-phase-based approach to 
macroscopic polarization of dielectrics has come to be known as the "modern theory of po- 
larization" (MTP) in the current physics literature .0 The MTP has been used both within 
the density-functional theory (DFT) based implement at ions, !'Ul!0~i! as well as the wave- 
function-based Hartee-Fock (HF)ir@ formulations, to evaluate a variety of polarization 
related properties. 

DFT-based calculations of macroscopic polarization properties are very efficient, so that, 
without excessive effort, one can perform ab initio computations on complex compounds. 
Normally, the ab initio Born charges computed using the DFT-based formulations are within 
10% agreement with the experiments for simple zinc-blende semiconductors^ while the dis- 
agreement can be worse for more complex systems.0 Therefore, it is of interest to system- 
atically explore alternative methods for computing the macroscopic polarization properties 
of crystalline insulators. For ionic systems, the HF method provides a powerful alternative 
in that it performs much better than the local-density approximation based schemes.iIII 
The other advantage of the HF method is that it can be systematically improved, both 
by perturbative, as well as nonperturbative methods, to account for many-body effect s.0 
Recently, we have developed a Wannier-function-based ab initio HF approach to compute 
the ground-state properties of crystalline insulators.iH§ The approach has been successfully 
applied to compute the ground-state properties of not only three-dimensional crystals,EjE3 
but also of one-dimensional periodic insulators such as polymers.lHiiil In the present paper 
we extend our approach to the problem of macroscopic polarization of dielectrics, and apply 
it to obtain benchmark HF values for the Born charges of diatomic ionic systems LiH, LiF, 
LiCl, NaF and NaCl. Besides their simplicity, the main criteria behind the choice of these 
materials for the present study were: (a) to the best of our knowledge, no prior benchmark 
calculations of the Born charges of these materials exist (b) high-quality experimental data 
has been available for them for a long time.llrll Thus by comparing the benchmark HF 
results such as this one, to the corresponding experimental values, one can gauge the appli- 
cability of the HF method on a wide variety of systems. When we compare the HF values 
of the Born charges computed in the present work, with the experimental ones, we find that 



2 



the agreement is good for LiH, LiF, and NaF, where the agreement between the theory and 
the experiment is always better than 7%. For NaCl and LiCl, however, the error is 10% and 
16%, respectively, suggesting the possibility that the many-body effects may be stronger in 
these systems. 

Since, this is the first application of our Wannier-function-based method to the problem 
of macroscopic polarization, we also present the associated computational details. The 
Wannier functions, being very similar in character to the molecular orbitals encountered 
routinely in quantum chemistry, have the added advantage of being intuitive in character. 
Indeed, as demonstrated later on, they lead to a pictorial description of the symmetry 
breaking processes associated with macroscopic polarization of insulators. 

Remainder of the paper is organized as follows. In section |T| we briefly cover the theo- 
retical aspects of the present work. Our numerical results for the Born effective charges of 
several ionic crystals are presented in section |HI[ Finally, the conclusions are presented in 
section [IV]. 



II. THEORETICAL BACKGROUND 

A. Born Effective Charge 

The Born effective charge tensor, Z*Ji), associated with the atoms of the i-th sublattice, 
is defined as0 



Qp(el) 



Z* a p(i) = Z i+ (n/e,— 



i/3 



(1) 



E=0 



where Zi is the charge associated with the nuclei (or the core) of the sublattice, Q is the 
volume of the unit cell, e is the electronic charge, and Pjf 1 ^ is the a-th Cartesian component of 
the electronic part of the macroscopic polarization induced as a result of the displacement of 

r\ p(eZ) 

the sublattice in the /3-th Cartesian direction, mp. For small Au^, one assumes 



i/3 



E=0 



A ^ } , and computes the change in the polarization APjf^ following Resta's approach^ 

AP el = P« - Pg> , (2) 

where, and P *p , respectively, denote the electronic parts of the macroscopic polarization 
of the system for its initial (A = 0) and final (A = 1) states, where A is a parameter 
characterizing the adiabatic symmetry-breaking transformation of the lattice. Clearly, for 
the present case, A is to be identified with the sublattice displacement A-u i/3 . For one-electron 
theories such as the Kohn-Sham theory, or the HF theory, King-Smith and Vanderbilt showed 
that@ 

M r 

P^ = (/e/fi)E/r4 A) (r) 2 rfr, (3) 

n=l J 

where { (p^(r), n — 1, . . . , M } represent the M occupied Wannier functions of the unit cell 
for a given value of A, and / is the occupation number of Wannier function (/ = 2, for the 
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restricted-Hartree-Fock theory). King-Smith and VanderbiltlI3 also showed that the r.h.s. of 
Eq. ([|) is proportional to the sum of the Berry phases associated with individual Wannier 
functions (or bands), thus equating the change in the macroscopic polarization of the solid 
with the change in the Berry-phase of its wave function during the corresponding adiabatic 
transformation. In addition, Resta0 demonstrated that AP e ; computed via Eqs. (0) and 
(0) is invariant under the choice of Wannier functions, even though the individual P^'s are 
not. Computation of the Wannier functions for different values of A is discussed in the next 
section. 



B. Wannier Functions 

In principle, any approach which can yield Wannier fucntions of a crystal corresponding 
to its Bloch orbitals, can be used to compute its Born charge tensor. However, in the 
present work we have applied a framework, recently developed by us, which directly yields 
the restricted-Hartree-Fock (RHF) Wannier functions of a crystalline insulator employing an 
LCAO approach.liE In our previous work we showed that one can obtain M RHF Wannier 
functions, {\oi), a = 1, Mioccupied by 2M electrons localized in the reference unit cell C 
by solving the equationsilll 

(T + U + J2(2Jp + A*| 7 (H*)><7(H*)I)I«> = (4) 

P keAf 7 

where T represents the kinetic-energy operator, U represents the interaction of the electrons 
of C with the nuclei of the whole of the crystal, while Jp, Kp, respectively, represent the 
Coulomb and exchange interactions felt by the electrons occupying the (3-th Wannier function 
of C, due to the rest of the electrons of the infinite system. The first three terms of Eq.(|J) 
constitute the canonical Hartree-Fock operator, while the last term is a projection operator 
which makes the orbitals localized in C orthogonal to those localized in the unit cells in 
the immediate neighborhood of C by means of infinitely high shift parameters A^'s. These 
neighborhood unit cells, whose origins are labeled by lattice vectors are collectively 
referred to as N ' . The projection operators along with the shift parameters play the role 
of a localizing potential in the Fock matrix, and once self-consistency has been achieved, 
the occupied eigenvectors of Eq.(^j) are localized in C, and are orthogonal to the orbitals of 
Af — thus making them Wannier functional!. As far as the orthogonality of the orbitals 
of C to those contained in unit cells beyond M is concerned, it should be automatic for 
systems with a band gap once M has been chosen to be large enough. As in our previous 
calculations performed on three-dimensional ionic insulators,"^ we included up to third- 
nearest-neighbor unit cells in the region N ' . 

For computing the Born charges, first Pg° is computed from Eq. (|3|) using the Wannier 
functions of the unit cell, with all the sublattices of the crystal in their original position 
corresponding to the case A = 0. Next, the z-th sublattice is displaced in the Cartesian 
direction (3 by a small amount Auip, and P^ is computed in a manner identical to the 
previous case, except that the Wannier functions used for the purpose are recomputed for 
the transformed lattice. Now that we can compute AP e ; (cf. Eq. (fj)), and compute the 
Born effective charge tensor by substituting it in Eq. ([!]). 
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The Wannier functions obtained by solving Eq. (Q) are canonical Hartree-Fock solutions 
for the unit cell C, and thus will satisfy the spatial symmetries of the unit cell. In appearance 
they look identical to the molecular orbitals encountered in any quantum-chemical calcula- 
tion on a finite system, as was discussed in our earlier work.Hll Therefore, by comparing the 
spatial appearances of the Wannier functions for the most symmetric case (A = 0), to the 
broken symmetry one (A = 1), we can obtain a pictorial representation of the polarization 
process. 



III. RESULTS AND DISCUSSION 

In this section we present the results of our calculations of the Born effective charges 
for LiH, LiF, LiCl, NaF, and NaCl. Because of the cubic nature of the underlying Bravais 
lattices, the Born charge tensor for these systems has only one independent component. In all 
the cases, we assumed the corresponding experimental fee crystal structure, with the anion 
at (0,0,0) position and the cation at (a/2,0,0) position, a being the lattice constant. The 
lattice constants used in the calculations were the theoretical ones, obtained by minimizing 
the total energy per unit cell at the Hartree-Fock level. 

The Wannier functions used in the approach were obtained by performing all-electron 
HF calculations using a computer program developed by us recently.El In order to evaluate 
the centers of the Wannier functions needed for computing the polarization properties, we 
added a small subroutine to the existing module. The program is implemented within an 
LCAO approach, employing Gaussian lobe-type functions.c2l Lobe-type functions simulate 
the Cartesian p and higher angular momentum orbitals located on a given atomic site, as 
linear combinations of s-type functions slightly displaced from the site.@ Because of this 
reason, it is possible, that in our approach we obtain somewhat different numerical values 
of the Wannier function centers, as compared to the ones computed by equivalent genuine 
Cartesian basis functions as implemented, e.g., in the CRYSTAL95 program^. For these 
calculations we used the lobe representation of the state-of-the-art contracted Gaussian basis 
sets developed by the Torino groupSS For LiH, the details of the basis set can be obtained 
in ref.e3, while for the alkali halides, they are available in ref.". Using these basis sets we 
had studied NaClH earlier at the Hartree-Fock level, therefore, optimized lattice constants 
were already available for them. However, for the remaining systems, we performed fresh 
Hartree-Fock calculations to obtain the optimized lattice constants. The theoretical lattice 
constant finally used in these calculations were 4.106 A (LiH), 4.018 A (LiF), 4.633 A (NaF), 
5.262 A (LiCl), and 5.785 A (NaCl). These are in close agreement with the values 4.102 A 
(LiH)§, 4.02 A (LiF)@, 4.63 A (NaF)@, 5.28 (LiCl)@, and 5.80 A (NaCl)@ reported by the 
Torino group. 

The computed Born effective charges are presented in table |. These results were obtained 
by translating the sublattices of a given crystal by the amount Au = 0.01a in the (100) 
direction. However, in order to ensure the stability of the results, several calculations were 
performed with different directions and magnitudes of Am, and no significant changes in 
the results were observed. It was also verified by explicit calculations that the sum total of 
all the effective charges corresponding to the different atoms of a unit cell was always zero, 
in agreement with its electrical neutrality. All these reasons give us confidence as to the 
correctness of our results. 
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From table | it is obvious that the theoretical Born effective charges obtained for LiH 
and the alkali halides are quite close to their nominal ionicities. This is in perfect agreement 
with the intuitive picture of these systems being highly ionic in nature. As far as the 
comparison of the HF results with the experimental results is concerned, it is very good 
for LiH, LiF, and NaF. However, for NaCl and LiCl, the disagreement is more than 10%. 
Similar differences with respect to the experiments were also observed by Yaschenko et al.Ea 
who computed the HF Born charge of MgO to be 1.808, while the experimental value for 
that compound is in the range 1.96 — 2.02.O One possible reason for the discrepancy between 
the theoretical and the experimental values of the Born charges could be the missing many- 
body effects. A qualitative discussion of these many-body effects was given by Harrison, in 
the context of his "ion-softening theory" .il When, e.g., the anionic sublattice of an alkali 
halide is translated, the bulk of the contribution to the Born charge — which, for the HF 
case, we call the mean-field contribution — is due to the electron transfer along the direction 
of the movement of the anion, and is associated with top-most occupied p-type Wannier 
function, as is obvious from Fig. [l|. However, according to Harrison^!, because of the many- 
body effects, we can have a single (virtual) excitation from the top p-type occupied Wannier 
function (the bonding orbital) into the first unoccupied Wannier function (the antibonding 
orbital) on the nearest-neighbor cations, thereby, modifying the Born charge. This virtual 
charge fluctuation, in effect, introduces some covalency into the system as compared to the 
mean-field HF results. In its simple parametrized form, the ion-softening theory of Harrison 
predicts a uniform value of Z* = 1.16, for the alkali halidesS This value of Z*, although 
reasonable, is clearly at variance with the experimental results which show a clear variation 
in the Z* values of different alkali halides. Therefore, it is of interest to borrow the essence 
of the many-body effects incorporated in the ion-softening theory, and apply it to these 
systems within a rigorous ab initio formalism, to test its applicability. Indeed, this is what 
we intend to explore in a future paper. 

In table |IJ we give the detailed contributions of various Wannier functions to the Born 
effective charges of the alkali halides, when the anionic sublattice is translated. It is clear 
from the table that the low-lying core-like orbitals basically translate rigidly along with the 
nuclei. Nonrigid translation is seen mainly for the ns and np Wannier functions of the anion, 
where n defines the top of the valence band. In particular, ns orbital gains some effective 
charge at the expense of the np orbital. The case of NaF is an exception to this rule where 
the Na 2p Wannier function makes a significant contribution to the effective charge (-0.216). 
However, this contribution is due to an accidental near degeneracy of the sodium 2p Wannier 
function with the 2s Wannier function of fluorine, which leads to their mixture when the 
HF equations (cf. Eq. (H)) are solved. Because of this reason, some of the Born effective 
charge associated with the 2s Wannier function of F is transferred to the 2p function of Na 
(cf. table H). However, this is an instructive example of the nonuniqueness of the individual 
Wannier functions. But, as should be the case, the total Born charge of fluorine in NaF 
is free of this ambiguity associated with the individual Wannier functions, in that it has a 
normal value of -0.956. 

It is also instructive to examine the polarization process pictorially, as depicted by Wan- 
nier functions. We shall do so for the specific case of NaCl. The Is and 3p Wannier functions, 
localized on the Cl~ site of the unit cell, are plotted along the (100) direction, in Figs. || 
and [I], respectively, both before, and after, the translation of the CI sublattice. As discussed 



6 



earlier, on intuitive grounds we would expect the highly localized Is Wannier function, which 
is the deepest-lying core orbital, to move rigidly with the nucleus. On the other extreme, 
we would expect the 3p Wannier function, which forms the top of the valence band, to show 
significant nonrigid behavior, because of its relatively delocalized character. This is indeed 
what we observe in Figs. and [I], respectively. Owing to the perfectly cubically-symmetric 
crystal field that the CI site sees in the undeformed lattice, one would expect the corre- 
sponding 3p Wannier function to exhibit perfect antisymmetry about its center. Once the 
CI sublattice is moved along the (100) direction, the crystal symmetry is reduced, and one 
would expect to see the signatures of the broken symmetry in the 3p Wannier function of CI. 
Both these phenomenon are clearly visible in Fig. [1], where, for the undeformed lattice the 
3p Wannier function is perfectly antisymmetric about its center, while for the deformed case, 
it is no longer so, and it shows clear signs of induced polarization due to broken symmetry. 



IV. CONCLUSIONS 



In conclusion, we have applied the Berry-phase-based theory of macroscopic polarization, 
developed by King-Smith and VanderbiKo, to obtain the benchmark values for the Born 
effective charges of several ionic compounds at the Hartree-Fock level. In the present work, 
we have utilized the Wannier functions as the single particle orbitals, and demonstrated 
that they lead to a pictorial description of the polarization process. As far as our results 
are concerned, they are in good agreement with the experiments for all the systems except 
LiCl and NaCl, where the disagreement with the experiments was more than 10%. One of 
the reasons behind this disagreement could be that the many-body effects in these systems 
are significant. Although, there have been generalizations of the theory of macroscopic 
polarization to include many-body effects,0 their implementation is not as straightforward 
as the single particle theory. Recently, we have generalized our Wannier-function-based 
approach to include many-body effects by systematically enlarging the many-particle ground- 
state wave function by considering virtual excitations from the space of the occupied Wannier 
functions to that of the virtual ones.il The approach was demonstrated by computing the 
correlation contributions to the total energy per unit cell of bulk LiHUH. In a future paper, 
we intend to generalize our approach to compute the influence of many-body effects on 
macroscopic polarization properties as well. 
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FIGURES 



1 




-1 1 1 1 1 1 1 1 

-6 -4 -2 2 4 6 

r (atomic units) 

FIG. 1. Wannier functions corresponding to one of the 3p valence orbitals of CI - in NaCl before 
(solid line) and after (dashed line) the CI sublattice translation, plotted along the (100) direction. 
The deformed lattice was obtained by translating the CI sublattice by 0.01a in (100) direction. 
Unlike the core orbitals (see Fig. ^), the valence Wannier function translates with significant 
nonrigid character, and shows signatures of broken symmetry. 




r (atomic units) 



FIG. 2. Wannier functions corresponding to the Is core orbital of Cl~ in NaCl for the unde- 
formed lattice (solid line), and the deformed lattice (dashed line). Rest of the information is same 
as in the caption of Fig. ||. Clearly, as expected, the Is Wannier function translates rigidly with 
the sublattice. 
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TABLES 



TABLE I. Born effective charges of different ionic crystals obtained in this work, as compared 
to the experimental values. 



Crystal 


Born Effective 
This Work 
Z* (theory) 


Charge 

Experiment 

Z*[exp) 


Percentage Error 

{Z* {theory)- Z*{exp)) „ 10 „ 
Z* (exp) 


LiH 


1.046 


0.991±0.04 a 


5.5 


LiF 


0.998 


1.045 fe 


-4.5 


LiCl 


1.036 


1.231 fe 


-15.6 


NaF 


0.956 


1.024 fe 


-6.6 


NaCl 


0.985 


1.099 6 


-10.4 



a Obtained from experimental values of the Szigeti charge Z s , and the high-frequency dielec- 
tric constant e^, using the relation Z* = ( £ °°^+2) Z S M Experimental Z s and were reported 
in Ref.ll 

b Obtained from the experimental values of Z s and reported in ref.il 



TABLE II. The contribution of individual Wannier functions to the Born effective charge of 
the alkali halides when the anion (A) sublattice with nuclear charge Z nuc was translated, holding 
the cation (C) sublattice(s) fixed. Nominal ionicity of the anion is given in the parenthesis right 
below its Born effective charge. 



Wannier function 


Nominal Charge 




Born Effective 


Charge 








LiF 


LiCl 


NaF 


NaCl 


Is (C) 


0.000 


0.001 


0.002 


0.000 


0.000 


2s (C) 


0.000 






-0.007 


-0.001 


2 P (C) 


0.000 






-0.216 


0.002 


Is (A) 


-2.000 


-2.000 


-2.000 


-2.000 


-2.000 


2s (A) 


-2.000 


-2.061 


-2.002 


-1.830 


-2.001 


2 P (A) 


-6.000 


-5.938 


-5.998 


-5.903 


-5.999 


3s (A) 


-2.000 




-2.107 




-2.081 


3 P (A) 


-6.000 




-5.931 




-5.905 


Znuc 




9.000 


17.000 


9.000 


17.000 


Total 




-0.998 


-1.036 


-0.956 


-0.985 






(-1.000) 


(-1.000) 


(-1.000) 


(-1.000) 
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